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We present three methods for calculating the Feynman propagator for the non- 
relativistic harmonic oscillator. The first method was employed by Schwinger a half 
a century ago, but has rarely been used in non-relativistic problems since. Also 
discussed is an algebraic method and a path integral method so that the reader can 
compare the advantages and disadvantages of each method. 



The purpose of this paper is to discuss three methods for calculating the Feynman propa- 
gator. The methods are applied to the harmonic oscillator so that they can be compared by 
the reader. The harmonic oscillator was chosen because of its intrinsic interest and because 
it is the simplest non-trivial system after the free particle (see for instance references QQ) . 

The first method we will discuss was developed by Schwinger 3 to treat effective actions in 
quantum electrodynamics and is based on the solution of the Heisenberg operator equations 
of motion. The use of proper operator ordering and the subsidiary and initial conditions 
yields the propagator. The second method is based on algebraic techniques based on fac- 
torizing the time evolution operator using the Baker-Campbell-Hausdorff formulae.^ By 
using factorization, the completeness relations, and the value of the matrix element (x\p), 
we can determine the propagator. This method is close to the one presented in Refs. |6| and 
Q, but here we will use the Baker-Campbell-Hausdorff formulae in a slightly different way. 
The third method is a path integral calculation that is based on a recurrence relation for 
the product of infinitesimal propagators. As far as we know, this recurrence relation has not 
appeared in previous discussions of the one-dimensional harmonic oscillator path integral, 
although it is inspired by a similar relation in the three-dimensional system.— 

To establish our notation, we write the Feynman propagator for a time independent 
non-relativistic system with Hamiltonian operator H in the form: 



I. 



INTRODUCTION 
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where U(t) is the time evolution operator: 

U(t) = exp(-iHr/U), (2) 
and 9(t) is the step function defined by: 

{1 if r > 
(3) 
if r < 

II. SCHWINGER'S METHOD 

This method was introduced in 1951 by Schwinger in the context of relativistic quantum 
field theory 3 - and it has since been employed mainly in relativistic problems, such as the 
calculation of bosonio^ and fermioni o 10 ! 11 ' 12 ^ 3 Green's functions in external fields. However, 
this powerful method is also well suited to non-relativistic problems, although it has rarely 
been used in the calculation of non-relativistic Feynman propagators. Recently, this subject 
has been discussed in an elegant way by using the quantum action principle.— Before this 
time, only a few papers had used Schwinger's method in this context . 1 ^ 1 ^ 1 ? We adopt here 
a simpler approach that we think is better suited for students and teachers. 

First, observe that for r > 0, Eq. leads to the differential equation for the Feynman 
propagator: 

ih—K(x", a/; r) = (x"\H exp(-^Hr) \x) . (4) 
or v n ' 

By using the general relation between operators in the Heisenberg and Schrodinger pictures, 



d H (t) = e ^t/hQ se -iHt/n ^ (5) 

it is not difficult to show that if \x) is an eigenvector of the operator X with eigenvalue x, 
then it is also true that 

X(t)\x,t) = x\x,t), (6) 

where 

X(t) = e iHt/h^ e -iHt/h ^ ( ? ) 

and \x,t) is defined as 

\x,t) = e i&t/h \x). (8) 
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Using this notation, the Feynman propagator can be written as: 

K(x", x'; t) = (x", t\x', 0) , (9) 

where 

X(r)\x",r) = x"\x",t) (10a) 
X(0)\x',0) = x'\x',0) (10b) 

The differential equation for the Feynman propagator, Eq. (J3J), takes the form 

ih^-(x",r\x',0) = (x",t\H\x',0) . (r>0) (11) 

The form of Eq. (fTTj) is very suggestive and is the starting point for the very elegant 
operator method introduced by Schwinger.— The main idea is to calculate the matrix element 
on the right-hand side of Eq. (fTTj) by writing H in terms of the operators X{r) and X(0), 
appropriately ordered. Schwinger's method can be summarized by the following steps: 

1. Solve the Heisenberg equations for the operators X(r) and -P(r), which are given by: 

ih^X{t) = [X(t),H] and ih—P(t) = [P(t),H] . (12) 
Equation (|12j) follows directly from Eq. (j5J). 

2. Use the solutions obtained in step (1) to rewrite the Hamiltonian operator H as a 
function of the operators X(0) and X{r) ordered in such a way that in each term of H, 
the operator X{r) must appear on the left side, while the operator X(0) must appear 
on the right side. This ordering can be done easily with the help of the commutator 
[X(0), X(r)] (see Eq. (|2*K|) ). We shall refer to the Hamiltonian operator written in this 
way as the ordered Hamiltonian operator H or ^(X (t) , X (0)) . After this ordering, the 
matrix element on the right-hand side of Eq. (fTTj) can be readily evaluated: 

(x",t\H\x',0) = (xV|#o r d(x(r),X(0)^|x',0> 

= H(x",x';t)(x",t\x',0) , (13) 

where we have defined the function H. The latter is a c-number and not an operator. 
If we substitute this result in Eq. (Ill}) and integrate over r, we obtain: 

(x", t\x', 0) = C(x", x') exp J H(x", x'; r')dT^j , (14) 

where C(x",x') is an arbitrary integration constant. 
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3. The last step is devoted to the calculation of C(x",x'). Its dependence on x" and x' 
can be determined by imposing the conditions: 

d 

(x",t\P(t)\x',0) = -ih—(x",T\x r ,0) (15a) 

ox" 

(x", r\P(0)\x', 0) = +ih^- (x", t\x', 0) . (15b) 

ox' 

These equations come from the definitions in Eq. (jlOj) together with the assumption 
that the usual commutation relations hold at any time: 

[X(r),P(r)] = [X(0),P(0)}=th. (16) 

After using Eq. (|15|) . there is still a multiplicative factor to be determined in C(x", x'). 
This can be done simply by imposing the propagator initial condition: 

lim (x", t\x', 0) = 8{x" - x') . (17) 

Now we are ready to apply this method to a large class of interesting problems. In 
particular, we shall calculate the Feynman propagator for the harmonic oscillator. 
The Hamiltonian operator for the harmonic oscillator can be written as 

H = + -muj 2 X 2 (r) , (18) 

or 

H = ^l + - mu j 2 X 2 (0) J (19) 

because the Hamiltonian operator is time independent, despite the fact that the operators 
.P(t) and X(r) are explicitly time dependent. It is matter of choice whether to work with 
the Hamiltonian operator given by Eq. (|18|) or by Eq. (J19)) . For simplicity, we choose the 
latter. 

As stated in step (1), we start by writing down the corresponding Heisenberg equations: 



s*w = ^ < 2te > 

^-P{t) = -muj 2 X(t), (20b) 



whose solutions permit us to write for t = t that: 

X(t) = X(0)coscur + ^-Uincur. (21) 

muj 

For later convenience, we also write the corresponding expression for -P(r): 

P(t) = -muX(0)smoJT + P(0)coscur. (22) 

To complete step (2) we need to rewrite -P(O) in terms of X(r) and X(Q), which can be done 
directly from Eq. (jUJ: 



P(0) = — — - \X(t) - X(0) cos lot] . (23) 
sm[ujT) 

If we substitute this result into Eq. (JTHJ), we obtain: 



muj 2 



H = \X 2 {t) + X 2 (0) cosVr) - X(0)X(t) cos(cur) 

2 sin (cur) 

- X(t)X(0) cos(cur)] + ^mu 2 X 2 (0) . (24) 



Note that the third term in Eq. (j24j) is not written in the appropriate order. By using the 
commutation relation 

\X(0),X(t)] = \X(0),X(0)cos(ut) + -^-sm( ujt)] 

mui 

" l sin(cur), (25) 



mid 

it follows immediately that 

X(0)X(t) = X(r)X(0) + — sin cut. (26) 

mco 

If we substitute Eq. into Eq. (J22J), we obtain the ordered Hamiltonian: 

H ord = ~ . T( : [X 2 (r) + X 2 (0) - 2X(r)X(0) cos(cur)] - — cot(cur) . (27) 
z sin (^r J * 

Once the Hamiltonian operator is appropriately ordered, we can find the function 
H(x", x'; r) directly from its definition, given by Eq. (fT3*j) : 

H(x",x',t) 



(x", r| 


H\ 


x',0) 


(x", 


t\ 


x',0) 



— — [(x" 2 + x' 2 ) csc 2 (cur) — 2x"x cot(a;r) csc(u;r)] — — — cot (cur) . (28) 

_ 
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By using Eq. (|14|). we can express the propagator in the form: 

(x",t\x',0) = C(x",x') exp|-^y r rfr , [^((x ,/2 + x ,2 )csc 2 ^r , 

— 2x"x' cot(cL'T ) csc cut') — — — cotc^r'] j>. (29) 

The integration over r' in Eq. (j2Uj) can be readily evaluated: 

(x ,r|x ,0) = exp<^ . — - (x + x ) cob{ut )—2xx\>, (30) 

^sm{ujT) lztism{LUT) ) 

where C(x",x') is an arbitrary integration constant to be determined according to step (3). 

The determination of C(x",x') is done with the aid of Eqs. (jl5j) and (|17j). However, 
we need to rewrite the operators -P(O) and P(t) in terms of the operators X{r) and X(0), 
appropriately ordered. For -P(O) this task has already been done (see Eq. ()23j0 . and for -P(r) 
we find after substituting Eq. ^ into Eq. If22j) : 

P(r) = mu;cot(u;r) [X(r) - X(0) cos cur] - mcuX(0) sin(u;r) . (31) 

Then, by inserting Eqs. (|3*T|) and (|371|) into Eq. (jlfiaj) it is not difficult to show that: 

M^>=0. (32) 
dx" y 1 

Analogously, by substituting Eqs. and (j3Tf|) into Eq. ()15b|) we have that dC(x", x')/dx' = 
0. The last two relations tell us that C(x", x') = C, that is, it is a constant independent of x" 
and x' . In order to determine the value of C, we first take the limit r — > + on (x", t\x', 0). 
If we use Eq. (jHOJ), we find that: 



lim (x",t|x',0> = lim 4= exp^(x" - x') 2 ] = C\ —8(x" - x') . (33) 



If we compare this result with the initial condition, Eq. (jl7)l . we obtain C = y / mu/2iiih. 
By substituting this result for C into Eq. (J3*U|) . we obtain the desired Feynman propagator 
for the harmonic oscillator: 

K(x",x';t) = (x",r|x',0) 

■ / — 7 e*PS . — - (x +x cos wr -2ix . (34) 
2niri sin(wr) [ 2ft,sin(c<;r) L J J 
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In Sec. |V]we shall see how to extract from Eq. (|34|) the eigenf unctions and energy eigenvalues 
for the harmonic oscillator and also how to obtain, starting from the Feynman propagator, 
the corresponding partition function. 

For other applications of this method we suggest the following problems for the interested 
reader. Calculate the Feynman propagator using Schwinger's method for (i) the constant 
force problem; (ii) a charged spinless particle in a uniform magnetic field; and (iii) a charged 
spinless particle in a harmonic oscillator potential placed in a uniform magnetic field. 

We finish this section by mentioning that Schwinger's method can be applied to time- 
dependent Hamiltonians as well . 1 ^ 1 ^ It also provides a natural way of establishing the mid- 
point rule in the path integral formalism (see Sec. IIV)1 when electromagnetic fields are 
present.-^ 

III. ALGEBRAIC METHOD 

The origin of the algebraic method dates back to the beginning of quantum mechan- 
ics, with the matrix formulation of Jordan, Heisenberg, and Pauli among others. Here, we 
present an algebraic method for calculating Feynman propagators which involves manipu- 
lations of momentum and position operators . 4ll8ll9l2 ° l21 This is a powerful method because 
it is connected with the dynamical symmetry groups of the system at hand. A knowledge 
of the underlying Lie algebra can be used to calculate eigenvalues without explicit knowl- 
edge of the eigenfunctions.-i'iS*^ It can also be used to calculate propagators for a wide 
range of problems . 6 i 22 i 23 i 24 i 25 i 26 A coherent-state version of the algebraic method for different 
problems has been discussed also . 23 ^ 27 Because the use of these mathematical tools can be 
a bit cumbersome at first reading, we prefer to explore a simpler version of this method, 
which is close to that in Ref. 0. For this purpose, the calculation of the propagator for the 
one-dimensional harmonic oscillator is excellent. 

The Hamiltonian operator H for a non-relativistic system can usually be written as a sum 
of terms involving the operators P and X which do not commute. Hence, the factorization 
of the time evolution operator U(r) = exp(— irH/h) into a product of simpler exponential 
operators involves some algebra. This algebra deals basically with the commutation rela- 
tions among these non-commuting operators, and uses formulae generically known as Baker- 
Campbell-Hausdorff (see Eq. (jHSJ)). The use of those formulae is the essence of the algebraic 
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method, because it is easier to calculate the action of these simpler exponential operators 
on the states \x) or \p), than to calculate the action on these same states of the original time 
evolution operator. The algebraic method can be summarized by the following steps: 

1. First rewrite the evolution operator U (r) as a product of exponentials of the operators 
X, P, and PX. (Note that, in contrast with Schwinger's method, here the operators P 
and X are time independent, that is, they are in the Schrodinger representation.) The 
factorization can be done with the help of the Baker-Campbell-Hausdorff formula^ 

e A Be- A = C, (35) 

where A, B, and C are operators (for simplicity we omit the hat on the operators) 
and 

C = B + [A, B] + i[A [A, B}} + [A, [A, B]]] + . . . , (36) 
valid for any A and B. Equation ()36|) can be iterated as: 

C 2 = (e A Be- A )(e A Be- A ) = e A B 2 e~ A 

C n = e A B n e -A ^ 

If we expand exp(C) and identify each power C n in Eq. (|3T|). we find 

e c = eW^, (38) 

which can be inverted to give: 

e B = e~ A e c e A . (39) 

We then identify B = —irH/h and find a factorized form of the evolution operator for 
a conveniently chosen operator A. The specific choice for A depends on the explicit 
form of the Hamiltonian. This factorization can be repeated as many times as needed. 
Note that, in general, the operator C in Eq. (|36jh which is an infinite series with 
B and multiple commutators of A and B, is more complicated than the operator B 
alone, which is proportional to the Hamiltonian. However, if we choose the operator A 
conveniently, this series can terminate and the remaining terms from the commutators 
can cancel some of those terms originally present in B. A more systematic way of 
doing this factorization is to use the Lie algebra related to the problem under study. 
This way will be sketched at the end of this section. 
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2. Next substitute the factorized Hamiltonian into the definition of the Feynman propa- 
gator K(x", x'\ r), Eq. (JTJ), and calculate the action of the exponential of the operators 
X, P, and PX on the state \x). For the operator X this calculation is trivial and 
for P we just need to use the closure relation 1 = J dp|p)(p| and the matrix element 
(x\p) = (l/2irh) 1 / 2 exp(ixp). For the mixed operator PX we need to use the expres- 
sion: 

(p'\ expi-i-fPX/h) \p) = e~~<5(p' - e~» , (40) 

where 7 is an arbitrary parameter to be chosen later. Equation (jjOj) comes from the 
relation 

exp(-2 7 PX/ft)|p) = e- 7 |e-» . (41) 

Before we apply the algebraic method to a specific problem, we derive Eq. (|JTJ). We first 
note that 

exp(-i-fPX/h)P\p) = pexp(-i^PX/h)\p) . (42) 

Equation (|4*2j) can be rewritten as 

[exp(-i-fPx/h)Pexp(i-fPX/h)] exp(-ijPX /fi)\p} = pexp(-ijPX/h)\p) , (43) 

so that we can use the Baker-Campbell-Hausdorff formulae (|33j) to rewrite the term in the 
square brackets as: 

exp(-ijPX/h)Pexp(i-fPx/h) = (l + 7 + ^ 7 2 + ^ 7 3 + . . )P 

= e 7 P. (44) 

If we substitute the above result into Eq. (J43j) . we have: 

P[exp(-i 7 Pi:/^)|p)] = e- 7 p[exp(-i 7 PX/^)|p)], (45) 

which shows that exp(— i'yPX /h)\p) is an eigenstate of the operator P with eigenvalue pe -7 . 
This eigenstate can be written as |e~ 7 p), up to a constant C 7 , so that 

exp(-i>yPX /%)\p) = C^e^p) . (46) 

To determine the constant C 7 , we note that: 

(p'\ exp(^ 7 XP) exp(-ijPX /h)\p) = |C 7 | 2 (e"V|e-». (47) 



10 



If we use the relation [XP, PX] = 0, and Eq. f)H9|1 . we have: 



so that 



e^5(p-p') = |C 7 | 2 e 7 5(p'-p). 



(48) 



(49) 



Equation (fUJj) determines C 7 = e -7 and finally yields Eq. (}4*Tj) . 

Now we are ready to apply the algebraic method to solve some quantum mechanical 
problems. For the harmonic oscillator the time evolution operator (J2J) becomes: 



P 2 1 

U(r) = ex P [-ir(— + -mu 2 X 2 )/h] . 



(50) 



-2m 2 

We follow step (1), choose A = aX 2 where a is an arbitrary parameter and B = —irH/h, 
and obtain from Eqs. (jHSj) and ([3H]) 



exp(aX 2 ) exp(—irH /h) exp(— aX z 



-it 



exp 



P 2 iha.~-i - - . m 

— + (XP + PX) + - 

2m m 2 









l>-( 2ah Y" 




j 


\ m J 







(51) 



Note that even though the Baker-Campbell-Hausdorff formulae have an infinite number of 
terms, the number of non-vanishing commutators between H and X 2 is finite. 

The general idea of the algebraic method is that we want to factorize the time evolution 
operator. In this case, a step in this direction corresponds to canceling the term with X 2 in 
the right-hand side of Eq. (J51|) . This is easily achieved with the choice 

muj 



a 



2h 



(52) 



Then by using the commutation relation [X, P] = ih, we have: 
exp(-irH/h) = e luJT/2 exp(-aX 2 ) exp 



i P 2 

h V 2m ' 



exp(aX' 



(53) 



We can repeat step (1) to reduce the above operator containing P 2 between brackets into 
a product of simpler terms. This time we need to use products of P 2 instead of X 2 . If we 
use Eqs. fl33j) and (j3Hj). we obtain: 

P 2 „ I ... 

exp(/3P 2 )(— + iuoPX) exp(-l3P 2 ) = (— + 2ujh(5)P 2 + iuPX (54) 
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and to eliminate the term proportional to P 2 on the right-hand side we take: 



1 



(55) 



which gives from Eqs. (jHSj) and ()54|) 



exp[-^r( 



p2 



iujPX)~\ = exp {-(3P 2 ) exp[--(iur)PX] exp{(3P 2 



h 2m 

If we substitute Eq. into Eq. (J33J), we have that: 



exp( -^rtf ) = e^ r/2 exp(-aX 2 )exp(-/3P 2 )exp 



(56) 



--(iur)PX 



exp(/?P 2 ) exp(aX 2 ) . 



(57) 

Equation (JSTj) is the expression for the time evolution operator written as a product of 
simpler operators obtained by applying step (1) of the algebraic method. 

We next follow step (2), insert Eq. (J57j) into the definition of the Feynman propagator 
Eq. (JTJ), and find 

K [x , x ; T) = exp[— a[x —x ) H — —J 



>< , exp [~(p'x" - px') - f3{p' 2 - p 2 )} (p'\ exp [-Uiu)r)PX\ \p) . 



2nh 



(58) 



If we use Eq. f!40|) with 7 = iur, and the definitions ()52|) and ()55|) . we have 

1 



K(x",x';t) 



exp 



2nh 
x / dp exp 



Ilia.' I „ 2 ,2, . Je iWT x" -X'^ 2 

[x — x ) + 2- 



IUJT 



l-e 



-2iuiT 



l-e 



-2iuir 



Amuoh 



p — 2imuj- 



l-e 



-2io;r 



(59) 



This integral has a Gaussian form and can be easily done giving the harmonic oscillator 
propagator: 



K(x",x';t) 



muj 



27rihsmuT 



exp 



imuj 



2/i sin u;r 



/ n2 i2\ // / 

[X + X ) COS LOT — 2x X 



(60) 



where we used that (l — e~ 2luJT ) = 2ie~ luJT sin cur and Euler's formula, e luJT = cosur + i sin cur. 
This result naturally agrees with the one obtained in Sec. IIHI using Schwinger's method. 

Before we finish this section, we want to comment that the algebraic method can be 
discussed on more formal grounds, identifying the underlying Lie algebra, and using it to 
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explicitly solve the problem of interest. For the one-dimensional harmonic oscillator we can 
find a set of operators 

L- = -~d xx , L + = ^x 2 , L 3 = ^xd x + ^, (61) 

such that the Hamiltonian operator can be written as H = (K/m)L- + muj 2 L + . The above 
operators satisfy the SO (3) Lie algebra 

[L + ,L_] = 2L 3 , [L 3 ,L ± ] = ±L±. (62) 

This algebra is isomorphic to the usual SU(2) Lie algebra of the angular momenta and 
can be used to construct specific Baker- Campbell-Hausdorff formulae , 2 ^ 21 so that given the 
algebra, the solution arises naturally.^ If one considers three-dimensional problems, which 
are more involved because of the presence of terms proportional to 1/r 2 , this algebra can 
still be used to find the propagator, but the operators have to be generalized.—^ These 
generalized operators can be used to solve a wide range of problems . 22 i 23 i 24 i 25 i 26 



IV. PATH INTEGRAL METHOD 

The path integral formalism was introduced by Feynman^ in 1948, following earlier 
ideas developed by Dirac.-^ In the last few decades, path integral methods have become 
very popular, mainly in the context of quantum mechanics, statistical physics, and quantum 
field theory. Since the pioneering textbook of Feynman and Hibbs,^ many others have been 
written on this subject, not only in quantum mechanics,— i 31 i 32 i 33 i 34 but also in condensed 
matter,— as well as quantum field theory , 36 i 37 i 38 to mention just a few. 

In this section we shall apply Feynman's method to again obtain the harmonic oscillator 
quantum propagator already established in Sees. ITT1 and IIHI The purpose here is to evaluate 
the corresponding path integral explicitly, without making use of the semiclassical approach 
which is often used in the literature. Of course this kind of direct calculation already exists in 



the literature, see for example Refs. , |38|, and 39. However, we shall present an alternative 
and very simple procedure. 

The path integral expression for the quantum propagator is formally given by: 

K(x N ,x ;r) = [ [Dx)e iS W h , (63) 

x(t)=x n 
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where S(x) is the action functional: 



*JV 1 



S{x) = I [-mx 2 (t) - V(x(t))]dt, 



(64) 



to 



and [Dx] is the functional measure. According to Feynman's prescription, we have that: 

N-l 



K(x n ,x ;t) 



lim 

N—>oo 



x exp< 



rn 



2-nihe 



n 



m 

2mfie Xj 



N 

£ 

fe=i 



2e 



(65) 



where Ne = r. With this prescription, the scenario is the following: summation over all the 
functions x means to sum over all the polynomials in the plane (t,x(t)), starting at (xo,t ) 
and finishing at (xn^n), which gives rise to the integrations over the variables Xj = x(tj) 
from — oo to oo, where tj = t Q + je, with j = 1, 2, . . . , iV — 1. Hence, to evaluate a path 
integral means to calculate an infinite number of ordinary integrals, which requires some 
kind of recurrence relation. 

When electromagnetic potentials are absent as is the case here, it is not necessary to 
adopt the midpoint rule for the potential V(x) as given by Eq. (J65|) . and other choices can 
also be made. Instead of using the midpoint rule we shall write the discretized version of 
the action as 

S = Y1 T,±(V{.r,) + V{ .r , _ , ) ) . (()()) 



j= i J 



where for generality we have taken Tj as the jth time interval so that r = ^ Tj. Then, we 
write the Feynman propagator in the form: 



. N N-l 

K(xn, xq] t) = Jim / Y\ K(xj, Xj-i, Tj) J^J dx 



(67) 



j=l k=l 

where the propagator for an infinitesimal time interval is given by: 



K[Xj, Xj—i] Tj) 



2nifiTj 



cxp 



m[Xj — Xj-i) 
2^ 



-rj^[V(xj) + V(xj^ 



(68) 



If we use the Lagrangian for the harmonic oscillator, namely, 

L(x, x) 



-mx muo x 

2 2 



(69) 
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the infinitesimal propagator, given by Eq. (Joaj) . takes the form: 

imuj 1 



2ix7i \ cur, 



cxp 



2Tl UJTj 



2x 3 x 3 A (70) 



To calculate this infinitesimal propagator we now define new variables <pj such that: 



sin (pj = UJTj, 



(71) 



which implies that <pj = ujTj and cos<pj = 1 — ojtj/2. In fact, other variable transformations 
could be tried, but this is the simplest one that we were able to find that allows an easy 
iteration through a convolution-like formula. It is also helpful to introduce a function F: 

imuj 1 



1 



2-niTi V sin 



cxp 



■ (cos <f)(r] 2 + rj ) — 2rj r\ ') 



(72) 



2Ti sin i 

Then, using Eqs. (fTUj) . (J7TJ), and (fT2~j) . we can rewrite the harmonic oscillator propagator 
(|H7|) as: 

„ N N-l 

K(x N ,x ;r) = lim / TTf^-, Xj~ x \ 4>j) T\ dx k . 

17 j=l k=l 

The function F has an interesting property: 



(73) 



F( V ", rj; </>")F(t,, r,'; vVv = Hv", V; 0" + 0')- (74) 
This can be seen by a simple direct calculation. ^From its definition (fT2*|) . we have that: 

/ F(v",V'A") F (v,v''A')d , r] = — IJTZIZ A / ex P 



2nih V sin 6" sin 



imui f cos 0" //2 cos </>' /2 



2fr V sin 



-77 + 



sin 



x / <ir7 exp 



2h 



■{arj 2 - 277/3) 



(75) 



where we have defined: 



sin(0" + 0') 
sin 0" sin <b' 



and /3 



77" sin <fi' + 77 ' sin 0" 



(76) 



sin 0" sin 0' 

By completing the square in the integrand of Eq. (J7HJ) and calculating the Fresnel integral, 
we obtain: 



F( V ", m 4>")F( V , V '; 4>')dv 



muj / 2-Kih 



2nih V mua V sin 0" sin 
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x exp 



imui ( COS 6 „2 COS0 , 2 3 

— r~ V H V 

2n V sin d>" sin ra' a 



(77) 



If we use the definitions in Eq. (176(1 . as well as some trivial manipulations with trigonometric 
functions, it is straightforward to show that: 



imuj 



-(cos{<j>" + 0')(V' 2 + V 2 ) - 2rf' V ')] ,(78) 



which is precisely Eq. (174(1 . We now define x" — x^, x' — xq, use the fact that 

lim^ 4>j — UJT an d the result in Eq. (|75j) to finally obtain the desired propagator: 



K(x",x';t) = F(x",x';ujt) 



2nih sintur 



exp 



imuj 



( ill i /2n // / 

cosu;t(£ +x ) — ix x 



(79) 



2% sin cur 

This result coincides with the expressions obtained in the previous sections for the harmonic 
oscillator. Let us now review some of the applications of the Feynman propagator. 



V. EIGENFUNCTIONS, EIGENVALUES AND THE PARTITION FUNCTION 

In this section we will show how to obtain the stationary states and the corresponding 
energy levels, as well as the partition function of the quantum harmonic oscillator, directly 
from the expression of the Feynman propagator calculated earlier. Although these tasks are 



38 and 



401 ) . we shall present them here 



well known in the literature (see for instance Refs. 
for completeness. 

To obtain the energy eigenstates and eigenvalues, we need to recast the propagator (|79|l 
in a form that permits a direct comparison with the spectral representation for the Feynman 
propagator given by: 

K(x, x'; r) = 6(r)]T ^{xW^e^ /h . (r > 0) (80) 

n 

If we define the variable z = e~ luJT , we can write: 

1 1 - z 2 

sin(a;r) = (81a) 

2% z 

l + z 2 

cos(cjt) = — — . (81b) 
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We further define £' = \J mu / Tlx' and £" = a/ muj/fix", and express the harmonic oscillator 
propagator to the form: 

K(x",x';r) = ^(1 - z 2 )^ 2 exp - + t" 2 ) (^f)] } (82) 



— (!-*) / e xp[--(e + £ )Jexp[ — J 



(83) 



where we used the identity 



2(1 - ^ 2 ) 2 I - : 
Now we consider Mehler's formula?^ 

~2xyz — (x 2 + y 2 )z 2 



1 + z 2 1 z 2 

~ +T^- (84) 



:i - ^r 1/2 



exp 



1-z 2 



Y^H n {x)H n {y)—. (\z\<l) (85) 



2 n n 

71=0 



However, some care must be taken in order to use Eq. (|83|) in Eq. (|82jl. because |z| = 1 and 
Mehler's formula (|85j) requires that |z| < 1. This problem can be circumvented if we add an 
imaginary part to u, namely, if we let u — > ui — z£, and take e — > after the calculations. 
Hence, if we use Eq. (J8*5)). Eq. (J52j) takes the form: 

^y iT) = ^ exp [ { ^ + ^ )] J : h^- X ")H^— X ')^^, 

V n=0 v v 

(86) 

where we have let x = £" = ^/rmj/hx" and y — ^' — a/ mu/Tix' . 

If we compare Eq. (jS^jl with the spectral representation ()80|). we finally obtain the well 
known results for the energy eigenfunctions (apart from a phase factor) and energy levels, 
respectively: 

<P n (x) = —=(—) 1 exp x *)H n L—x) 87 

E n = (n + ^)hu. (88) 

We finish this section by calculating the partition function for the harmonic oscillator. 
With this purpose in mind, recall that the partition function in general can be written as: 

Z(J3) = Tre- fi6 . (89) 
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The trace operation can be taken over a discrete basis, the eigenfunctions of the Hamiltonian 
itself, or, more conveniently here, over the continuous set of eigenstates of the position 
operator (denoted by \x)): 

/+oo 
dx(x\e- pH \x). (90) 
-oo 

If we identify the integrand with the Feynman propagator with the end points x' = x" = x 
and (3 = ihr as the imaginary time interval, we have 

/+oo 
dxK(x,x;-ihP). (91) 
-oo 

Then, from the harmonic oscillator propagator Eq. ()79j1 . we readily obtain: 



K(x, x; -ih/3) = J m< ? exp [- - , — (cosh(u /3h) - 1) x 2 ] , (92) 

where we have used sin(— ia) = — isinha and cos(ia) = cosh a. By substituting Eq. ([9*2^1 
into Eq. ()91|) . and evaluating the remaining Gaussian integral, we finally obtain: 

mux 2 ,uj@%y 



m = J tohJUiufih) 1 ( ' X|) 



oo 



oo 



tanh(— -) 



k v 2 



dx 



2sinh(i^) ' ^ 93 ^ 

where we used the identities cosh a — 1 = 2sinh 2 (a/2) and sinh(a) = 2 sinh(a/2) cosh(a/2). 
Equation (|93jl is the partition function for the one- dimensional harmonic oscillator. 



VI. CONCLUSIONS 



We have rederived the one-dimensional harmonic oscillator propagator using three dif- 
ferent techniques. First we used a method developed by Schwinger that is usually used in 
quantum field theory, but that is also well suited for non-relativistic quantum mechanical 
problems although rarely used. We hope that our presentation of this method will help it 
become better known among physics teachers and students. Then we presented an algebraic 
method that deals with the factorization of the time evolution operator using the Baker- 
Campbell-Hausdorff formulae. We hope that our presentation will motivate teachers and 
students to learn more about such powerful methods, which are closely connected with the 
use of Lie algebras. It is worth mentioning that these methods can be applied not only in 
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non-relativistic quantum mechanics, but also in the context of relativistic theories. Finally, 
we presented a direct calculation of the Feynman path integral for the one-dimensional har- 
monic oscillator using a simple but very convenient recurrence relation for the composition 
of infinitesimal propagators. We hope that the presentation of these three methods together 
with the calculation of the Feynman propagator will help readers compare the advantages 
and difficulties of each of them. 
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